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ABSTRACT 

Context. The Taurus-Auriga star-forming complex hosts the only population of T Tauri stars in which an anticorrelation of X-ray 
activity and rotation period has been observed. 

Aims. We aim to explain the origin of the X-ray activity-rotation relation in Taurus- Auriga. We also aim to put the X-ray activity of 
these stars into the context of the activity of late-type main-sequence stars and T Tauri stars in the Orion Nebula Cluster. 
Methods. We have used XMM-Newton's European Photon Imaging Cameras to perform the most sensitive survey to date of X- 
ray emission (0.3-10 keV) from young stars in Taurus-Auriga. We investigated the dependences of X-ray activity measures - X-ray 
luminosity, L x , its ratio with the stellar luminosity, L x /L*, and the surface-averaged X-ray flux, F X s - on rotation period and compared 
them with predictions based solely on the observed dependence of L x on a star's L+ and whether it is accreting or not. We tested for 
differences in the distributions of L x /L* of fast and slow rotators, accretors and non-accretors, and compared the dependence of L x /L+ 
on the ratio of the rotation period and the convective turnover timescale, the Rossby number, with that of late-type main-sequence 
stars. 

Results. We found significant anticorrelations of L x and F X s with rotation period, but these could be explained by the typically higher 
stellar luminosity and effective temperature of fast-rotators in Taurus- Auriga and a near-linear dependence of Lx on L*. We found 
no evidence for a dependence of L x /L* on rotation period, but for accretors to have lower L\/L+ than non-accretors at all rotation 
periods. The Rossby numbers of accretors and non-accretors were found to be the same as those of late-type main-sequence stars 
showing saturated X-ray emission. 

Conclusions. Non-accreting T Tauri stars show X-ray activity entirely consistent with the saturated activity of fast-rotating late-type 
main-sequence stars. Accreting T Tauri stars show lower X-ray activity, but this cannot be attributed to their slower rotation. 

Key words, stars: pre-main sequence - stars: activity - stars: rotation - X-rays: stars - Open clusters and associations: Individual: 
Taurus-Auriga 



1. Introduction 

T Tauri stars are young stars contracting toward the main se- 
quence, a subclass of which, the 'classical' T Tauri stars, show 
signatures that they are actively accreting material from circum- 
stellar accretion disks. Both accreting and non-accreting T Tauri 
stars exhibit strong, variable, X-ray emission which provides ev- 
idence for hot coronae generated by magnetic activity. Key open 
questions concerning this emission are the extent to which this 
magnetic activity is analogous to that exhibited by the Sun and 
similar stars, and what effect interaction with the circumstellar 
material, particularly through accretion, has on the magnetic ac- 
tivity. 

In the Sun, and all stars whose structure is composed of a 
radiative interior and a convective envelope, magnetic activity 
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is believed to be predominantly generated by a dynamo that is 
located in a shell at the interface of the radiative and convec- 
tive zones and is driven by (differenti al) rot ational and convec- 
tive motions (an aQ. dynamo: iParkerl Il955h . This is supported 
by observations of signatures of magnetic activity, such as chro- 
mospheric Ha or Ca H and K line emission or coronal X-ray 
emission, in stars with convective envelopes. The strengths of 
these activity signatures correlate directly with the projected 
rotation velocity, v s in ;', and inversely with the rotation period 
of the star, f rot (e.g. lPallavicini et al.L [l981). They furthermore 
show a tighter inverse correlation with the Rossby number, Ro, 
which is the ratio of the rotation period to the convective turnover 
timescale at the b ase of the convective envelope, T conv (e.g. 
iNoves etal.Lll984l) . 

The strengths of the activity signatures are observed to sat- 
urate at fast rotation velocities, short rotation periods and low 
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Rossby numbers dVilhul Il984t IVilhu & Waited 119871) . It is not 
yet known whether this is due to a saturation of the dynamo it- 
self or constraints on the available volume of the corona (see e.g. 
lJardine&UnruR[T999h . 

iPizzolatoe t al. (20031) have shown that the ratio of the X- 
ray luminosity to the stellar bolometric luminosity, Lx/L+, for 
main-sequence stars of spectral types from G to early M is 
well-determined by the Rossby number, Rq, such that Lx/£* ~ 
lO- 3 - 2 (Ro/R ,satT 2 for R > # ,sat and L x /L* * lO" 32 for 
Ro < ^o.sat, where /?o,sat ~ 0.2 is the Rossby number at which 
saturation sets in. 

T Tauri stars with masses greater than approximately 0.3 M 
are expected to be initially fully-convective and to develop a 
solar-like structure after nuclear fusion begins in the core. Stars 
with masses less than approximately 0.3 M Q are expected to re- 
main fully-convective. The majority of T Tauri stars are thus 
expected to have long convective turnover timescales of 100- 
300 d. Despite the fact that an a£l dynamo cannot operate in 
fully-convective main-sequence stars, there is no empirical ev- 
idence for different characteristics in the magnetic activity of 
these stars: fast rotators show activity at the saturated level and 
only slow rotators show activity well below th e saturated level 
dPelfosse et al.Lll998tlMohantv & Basrill2003l) . 

The T Tauri stars in Taurus-Auriga stand out as the only 
population of pre-main sequence stars in which an anticorrela- 
tion of X - ray ac t ivity with rotation period has been observed 
dBouyieti 1 19901: iDamiani & Micelal 119951: INeuhauser et all 
1 19951: IStelzer & Neuhauseri 1200 ll) . The largest survey to date of 
the X-ray emission from T Tauri stars, the Chandra Ultradeep 
Orion Project (COUP) observation of almost 600 T Tauri stars in 
the Orion Nebul a Cluster (ONC), found no such anticorrelation 
dPreibisch et al.L 120051) . The Rossby numbers of all the COUP 
stars with measured rotation periods placed them in the saturated 
regime of main-sequence stars, and their L x /^* was largely con- 
sistent with the saturation level of main-sequence stars but with 
a very large scatter. A strong dependence of X-ray luminosity 
on the stellar bolometric luminosity was found, as had been 
seen in previous surveys of X-ray emission from T Tauri stars 
(e.g. lPreibisch & Zinneckerll2002l) . Accreting T Tauri stars were 
found to have generally lower X-ray activity levels than non- 
accretors. 

The lower X-ray activity of ac creting T Tauri stars has 
also been observ e d in T a urus- Auriga (Damiani & Mi ceEll995t 
INeuhauser et all 119951: IStelzer & Neuhauseri |2001|) . and has 
been usually attributed to their typically slower rotation and an 
anticorrelation of X-ray activity with rotation period. 

We use the results of a new X-ray survey of Taurus-Auriga, 
the XMM-Newton Extended Survey of the Taurus Molecular 
Cloud (XEST), described in Sect. [2] to reinvestigate the depen- 
dence of X-ray emission on rotation in this star-forming region. 
Our sample selection and data analysis methods are described 
in Sects. [3] and |4] respectively. In Sect. [5] we present, and pro- 
pose an explanation for, the observed dependences of X-ray ac- 
tivity measures on rotation period. In Sect. [6] we explore the role 
of rotation in the different activity levels of accreting and non- 
accreting stars. We investigate the dependence of X-ray activity 
on Rossby number in Sect. [7] to examine the observed activity 
in the context of the activity of solar-like main-sequence stars 
and of T Tauri stars in the ONC studied by COUP. In Sect. |8] we 
make a careful comparison with, and reexamination of, the most 
comprehensive previous X-r ay study of Taurus-Auriga, made by 
IStelzer & Neuhauseri d200ll) with the ROSAT Position Sensitive 
Proportional Counter. In Sect. [9] we conclude by summarizing 



our findings and identifying outstanding questions and suggest- 
ing prospective observing strategies which could address them. 



2. The XMM-Newton Extended Survey of the Taurus 
Molecular Cloud (XEST) 

The XEST comprises 28 observations with exposure times of 
at least 30 ks of regions of the Taurus-Auriga complex that 
are most densely-populated by T Tauri stars. XMM-Newton' 's 
coaligned European Photon Imaging Camera (EPIC) PN, MOS 1 
and MOS2 detectors each have a full field of view approximately 
30 arcmin in diameter and perform simultaneous imaging, tim- 
ing and spectroscopy of X-rays in the energy range 0.2-12 keV, 
with respective resolutions of 6 arcsec (FWHM), less than 2.5 s, 
and 50-200 eV (E/AE m 20-50). The survey obser vations and 
data re duction procedures are described in detail in lGiidel et al.l 
d2007al) . 

XEST makes key improvements over previous surveys 
of Taurus-Auriga using the Einstein Imaging Proportional 
Camera (IPC) an d ROSA T Posi ti on Sensitive Prop o rtional 
Camera (PSPC ) dBouyierl 1 19901: [Damiani& Micelal 119951: 
INeuhauser et all 119951: IStelzer & Neuhauseri 1200 ll) thanks pri- 
marily to the greater sensitivity, broader energy band and higher 
spectral resolution of the EPIC detectors, augmented by longer 
exposure times. XEST reaches X-ray lumin osities of 10 28 erg s - 1 
at the 140 pc distance of Taurus-Auriga (Loinardet al., 2005), 
approximately an order of magnitude lower than t he ROSAT 
PSPC pointing survey (Stelzer & Neuhausetl 12001 1) , enabling 
the detection of almost all known T Tauri stars in the survey area 
down to the substellar mass limit for the first time. This is vital 
for accurately quantifying the dependence of X-ray activity on 
parameters such as stellar luminosity. T Tauri stars, particularly 
accreting T Tauri stars, are typically surrounded by circumstellar 
material or observed through molecular cloud material that ab- 
sorbs soft X-ray emission from the star. The broad bandpass of 
EPIC allows it to detect the less absorbed harder emission. The 
EPIC spectra enabled the absorption to be measured for each T 
Tauri star and accounted for in calculating its X-ray luminos- 
ity. This is crucial for accurate determination of X-ray luminosi- 
ties and was not generally possible with the data obtained by the 
Einstein and ROSAT surveys. Additionally, since the ROSAT sur- 
veys many more low-mass members of Taurus-Auriga have been 
identified by optical and near-infrared photometric surveys (e.g. 
iLuhman et all 120031) . making the known population of T Tauri 
stars much more complete. 



3. Sample selection 

We performed a literature study to compile a catalogue of recog- 
nised members of Taurus-Auriga and their relevant data, such 
as stellar luminosity, L*, effective temperature, T e s, multiplicity, 
photometric rotation period, P rot , projected equatorial rotation 
velocity, v sin i, and signatures of accretion such as the equiva- 
lent width of the Ha line. This catalogue and the many sources 
of information used in its compilati on are found, alon g with the 
XEST X-ray data for these stars, in lGiidel et all d2007al) . 

The XEST data provide an almost detection-complete sam- 
ple of more than 100 T Tauri stars (excluding Class I protostars, 
Herbig Ae stars and substellar objects) evenly divided between 
accretors and non-accretors and well-spread in masses from 0. 1 
to 2.5 M . We have chosen to restrict our study to T Tauri stars 
of spectral types M3 and earlier (r e ff > 3400 K). No star in the 
XEST sample of later spectral type has a measured photometric 
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period, and few such stars have measured v sin i. Stars of later 
spectral type are typically less-well studied and are not expected 
to develop a solar-like structure. 

We have excluded four undetected accreting T Tauri stars 
for which Lx upper limits could not be calculated due to lack 
of knowledge of absorption (Coku-Tau 1, HH 30, ITG 33 A and 
IRAS S04301+261), two T Tauri stars whose characteristic X- 
ray activity could not be assessed because a decay from a large 
flare persisted throughout their XEST observation (DH Tau and 
V830 Tau), and all T Tauri stars for which r e ff and/or L* were 
not available, or whose T e g and L+ placed t hem o utside the stel- 
lar evolutionary model grids of Si ess et ail d2000l) . The decay of 
a large flare also dominated the XEST observation of FS Tau, but 
we have used additional Chandra data to assess its characteristic 
X-ray activity. 

Our sample contains a total of 74 T Tauri stars, 47 accretors 
and 27 non-accretors. These include 23 stars with a measured 
photometric rotation period, of which 13 are accretors and 10 are 
non-accretors. There are 47 stars that have a measured vsirn, of 
which 30 are accretors and 17 are non-accretors, and this sample 
includes all stars with a measured photometric rotation period 
except for the accretor XZ Tau. 



4. Data Analysis 

Our study comprises four investigations. The data analysis per- 
formed in these investigations is described in this section, while 
the results of each investigation are presented and discussed in a 
dedicated section in Sects. [5}{8] 

Firstly, we examined the correlations with rotation period 
of the three most commonly-used measures of X-ray activ- 
ity: the X-ray luminosity, Lx, its ratio with the stellar bolo- 
metric luminosity, Lx/£*, and the surface-averaged X-ray flux, 
Fxs = Lx/4jtR+. In order to compare to previous studies we 
have treated accreting and non-accreting stars as a single sample. 
We tried to understand the observed activity-rotation relations in 
terms of the near-linear dependency of Lx on L+, and the typi- 
cally lower Lx of accreting stars, which have been consistently 
observed in populations of T Tauri stars tha t do not show an an- 
ticorrelation of activity on ro tation (e.g. iPreibisch et all 120051: 
Preib isch & Zinnecketi [2002). We used the correlations of Lx 
on L* derived s eparately for ac c retors and non-accretors in the 
XEST study by ITelleschi et ail d2007l see Fig. [B to calculate 
the expected Lx, and hence the Lx/L* and Fxs, of each star in 
our sample based on its and whether it was accreting or not. 
We compared the resulting correlations with rotation period of 
the expected and the observed activity measures. The results are 
presented and discussed in Sect. [5] 

Secondly, we tested the hypothesis that the lower Lx/L+ of 
accretors compared to non-accretors is due to their slower ro- 
tation period. In this analysis we have divided the full sample 
described in Sect. 3, into subsamples based on rotation period 
calculated from v sin i and on whether the star was accreting or 
not and statistically compared the distributions of Lx/L* of the 
subsamples. The results are presented and discussed in Sect. [6] 

Thirdly, we compared the X-ray activity of the T Tauri stars 
in our sample with the activity of solar-like main-sequence stars 
and T Tauri stars in the ONC based on the dependence of LxlL+ 
on Rossby number. The results are presented and discussed in 
SectH 

Finally we compared our activity-rotation relations with 
those obtained by the previous most-comprehensive X-ra y sur- 
vey of Taurus-Auriga, made bv lStelzer & Neuhauserl d2001l) us- 
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Fig. 1. The dependence of X-ray luminosity upon stellar lu- 
minosity for the accreting (black triangles) and non-accreting 
(grey circles) T Tauri stars in Taurus-Auriga observed by XMM- 
Newton. The black and grey lines show the best-fitting corre- 
lations for accretors and non-accr etors, respectively, fou nd us- 
ing the EM algorithm in ASURV ( ITelleschi et all 120071) . Open 
symbols and downward-pointing arrows mark upper limits. The 
cross shows KPNO-Tau 8, which was excluded from the fitting. 



ing the ROSAT PSPC. The results of this comparison are pre- 
sented and discussed in Sect. [8] 

The X-ray activity measures and rotation periods we used 
are described in Sects. |4~D and 14.21 respectively. The criterion 
we used to distinguish accretors and non-accretors is given in 
Sect. 14.31 The estimation of the convective turnover timescale 
for each of the stars in our sample is described in Sect. 14.41 The 
statistical methods used in our investigations are described in 
Sect.gT 



4. 1. Determination of the X-ray activity measures 

X-ray activity measures, Lx, Lx/L* a nd Fxs. were calcu lated 
for each star using the data tabulated in Gud el et all d2007al).The 
Lx of e ach observed T Tauri star was calculated in iGtidel etail 
d2007al) as follows. A spec tral fit of each source was made 
in XSPEC dArnaud , 1996) using a model composed of an 
optic ally- thin collisionally-ionized, so-called 'coronal', plasma 
(APEC: ISmith et all 200 1 ), having a two-power-law distribution 
of emission measure with temperature, multiplied with a photo- 
electric absorption model (WABS). The emission measure distri- 
bution was cut off at 10 6 and 10 s K. The elemental abundances 
of the plasma were fixed to values derived from high-resolution 
X-ray spectra of young starfl The power-law index at low ener- 
gies, a, is unconstrained in EPIC sp ectra and was fixed to + 2, as 
observed in young solar analogues dTelleschi et al.Ll2005h . The 
free parameters were the break temperature of the two power 
laws, 6.3 > log To > 7.5, the power-law index at high energies, 
-3 > fi > +1, the normalization, and the equivalent hydrogen 
column density of the absorbing gas, A^h- We integrated the flux 
under the best-fitting model between 0.3 and 10.0 keV and as- 



1 The abundance valu es used were, with respe ct to the solar pho- 
tospheric abundances of lAnders & Grevessd il989h : C=0.45, N=0.788, 
0=0.426, Ne=0.832, Mg=0.263, Al=0.5, Si=0.309, S=0.417, Ar=0.55, 
Ca=0.195, Fe=0.195, Ni=0.195. 



4 



K. R. Briggs: The X-ray activity-rotation relation in Taurus-Auriga 



sumed a distance of 140 pc to calculate the X-ray luminosity. 
Time intervals in which a star was clearly flaring were removed 
before spectral analysis if possible. 

We excluded observations of stars in which the emission was 
clearly dominated by a flare throughout the observation. If there 
were two or more acceptable observations of the same star we 
calculated the X-ray luminosity as the logarithmic average of 
the individual Lx determinations. 

For each observation of a star we estimated the uncertainty in 
Lx by finding the 68 per cent confidence interval in A^h and cal- 
culating Lx for the best-fitting models with A^h fixed to the lower 
and upper limits of this confidence interval. Lx was poorly con- 
strained when there were few counts, particularly when A^h was 
appreciable (> few xlO 21 ctrT 2 ) and the slope of the spectrum 
at energies above 2 ke V was not well-defined. A degeneracy be- 
tween log Tq and A^h sometimes developed, that in the worst case 
led to well-fitting models with very low log Tq and very high 
A^h, hence high Lx- Such cool temperatures were not observed 
in good-quality spectra showing low A^h and it is therefore likely 
that such models highly overestimate Lx- In such cases the er- 
ror bars plotted throughout this work extend from these probably 
unrealistically high values of Lx to lower values corresponding 
to more realistic hotter models with lower A^h- In some poorly 
constrained fits we fixed f3, usually th e least well-constr ained 
parameter, to a typical value of -1. In iGiidel et alj d2007al) we 
also fixed A^h in several particularly poorly-constrained cases 
to the value indicated by the visual or near-infrared extinction. 
However, this results in unrepresentative uncertainies and we 
have reanalysed those four cases here with just /3 fixed to - 1 . 
The results of this additional analysis are given in Table Q] 

We have also performed additional spectral analysis of 
four jet-driving accreting T Tauri stars, DG Tau A, GV Tau, 
DP Tau and CW Tau, whose X-ray spectra could not be well- 
described by a mo del with just a single absorbing column. 
IGiidel et all ((2007b) found instead a little-absorbed unvarying 
low-temperature component, which they proposed to be due to 
shocks in the jets, and a highly-absorbed and variable hot com- 
ponent, which they attributed to a corona. We have modelled 
the hot component with the spectral model describ ed above - 
rather th an the absorbed isothermal plasma used by Giide l et all 
(2007b) - and, the cool component with an isothermal model 
with photoelectric absorption. The hot components of DG Tau A 
and GV Tau were strongly variable in their XEST observations. 
In these cases we have first used the average spectrum from the 
whole observation to determine the best-fitting parameters of the 
unvarying cool component. The parameters of the cool compo- 
nent were then fixed to these values when we fitted the spectrum 
of the time intervals when the star appeared to be in a low or 
'quiescent' X-ray emission state to derive characteristic parame- 
ters for the hot component. CW Tau required no cool component 
when we excluded energies below 0.5 keV. We fixed /3 to -1 in 
these fits. Nevertheless, log Tq was unconstrained in all cases, so 
we have calculated Lx for log Tq = 7.0 and have expressed the 
range of uncertainties found when log Tq was a free parameter 
(see TableQ])- 

Just one star in our sample, FV Tau/c, was not detected in 
the XEST observation. FV Tau/c has no m easured photometri c 
rotation period or v sin i. As described in IGiidel et aU d2007al) . 
we calculated a 95 per cent confidence upper limit to its EPIC 
count-rate and used a spectrum with log Tq = 7.0 and f3 - -1, 
with A^h estimated from its optical extinction, to calculate an 
upper limit to its X-ray luminosity. 

For recognised multiple systems not spatially resolved into 
their individual components by XMM-Newton, we have assumed 



that the Lx of the components scales linearly with L*, as moti- 
vated by Fig. Q] Thus, each component of a multiple system is 
assumed to have the same Lx/L+, which is the ratio of the ob- 
served Lx and the total stellar luminosity of the system. The 
multiplicity-corrected Lx of the primary is therefore the ob- 
served Lx multiplied by the fraction of the summed stellar lu- 
minosities that is contributed by the primary. This fraction was 
estimated using the luminosities of the individual components if 
these had been determined, else the magnitude differences in the 
/f-band or /-band, as availably. Lx/£* was always calculated 
using the uncorrected observed Lx and the total stellar luminos- 
ity of the system. Fxs = Lxl^R\ = (Lx/ L^crT^ was always 
calculated using L x /L* and the effective temperature, T e ^ of the 
primary. 

4.2. Rotation periods 

Rotational data in the form of rotation per i ods an d spectroscopic 
v sin / have been taken from Rebull et al. (2004), where a com- 
prehensive list of references can be found, with a small num- 
ber of additions from L. Rebull (priv. comm.), and were always 
assigned to the primary in the cases of unresolved multiple sys- 
tems. For the larger sample of T Tauri stars with measured v sin i, 
we estimated rotation periods as P mt / sin i = 2nR*Jv sin i. These 
are essentially upper limits to the rotation perioqj but we note 
that for a randomly oriented sample (sin/) = 7r/4 » 0.785, the 
lcr lower limit to sin / is 0.73 and there is just a 10 per cent 
chance that sin i is less than 0.44. RY T au, LkCa 21 a nd CW Tau 
have tabulated photometric periods in Giid el et all (|2007a) but 
their high measured vsinZ indicate that t he photometric peri- 
ods a re too long to be the rotation periods (Bou vier et all 1993, 
1 19951) and we use only the v sin / values for these three stars in 
this work. 



4.3. Accretion criterion 

Accreting pre-main sequence stars, classical T Tauri stars, have 
been traditionally identified through their strong Ha emission, 
believed to arise from heated accreting material. Ha emission is 
also produced by chromo spheric activity, but, outside of excep- 
tional flares, is restricted by the saturation of magnetic activity 
and one can u se the criterion log(LH g /£*) > —3.3 t o identify 
accreting stars dBarrado v Navascues & Martini l2003h . It is sim- 
pler to measure the equivalent width of the Ha line, EWn a , but 
cooler stars have less continuum emission at the Ha line and 
so the criterion becomes a function of spectral type. We have 
adopted the dependenc e on spectral type of the EWn a criterio n 
that was proposed by Barra do v Navascues & Martin (20030, 
directly based on the criterion log(LH ff /L*) > -3.3. The clas- 
sification of the primary was used in cases of multiple systems 
not resolved by XMM-Newton. 



2 The X"-band is not ideal as accretion disks can emit substantially in 
this band, and lower-mass companions emit a higher proportion of their 
bolometric flux in this band relative to higher-mass primaries. For this 
latter reason, however, many multiplicity studies have been performed 
in this band. 

3 Two of the slowly-rotating non-accretors, HBC 358 and HBC 359, 
have only upper limits to v sin i of 10 km s and hence lower limits to 
P Iot / sin / of > 6.98 and > 6.72 d, respectively. 

4 The criterion proposed by Martfij i 19971) gives identical classifica- 
tions for the present sample. 
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XEST 


Name 


Nh (lo" range) 


logTo (lcr range) 


L x (lcr range) 


log(L x /L*) 




dof 






10 22 cm- 2 


K 


10 30 ergs" 1 






04-010 


GH Tau AB 


0.48 (0.33,0.63) 


6.4 (6.3,7.3) 


0.39 (0.29,0.50) 


-3.91 


0.89 


4 


09-010 


HO Tau AB 


0.27 (0.09,0.44) 


6.6 (6.3,6.9) 


0.064(0.034,0.13) 


-4.00 


1.11 


3 


07-011 


JH223 


0.19(0.14,0.31) 


6.6 (6.3,6.7) 


0.099 (0.078,0.18) 


-3.84 


0.89 


4 


11-079 


CFHT-Tau21 


1.72 (1.20,2.20) 


7.2 (6.3,7.5) 


0.22 (0.15,0.36) 


-3.82 


0.68 


7 


02-022 


DG Tau A 


4.50 (2.07,7.61) 


7.0 (6.3,7.5) 


0.55 (0.17,2.39) 


-4.08 


0.77 


13 


20-046 


CWTau 


2.93 (1.33,4.72) 


7.0 (6.3,7.5) 


0.11 (0.038,0.31) 


-4.60 


0.94 


9 


10-045 


DP Tau 


6.31 (2.98,9.02) 


7.0 (6.3,7.5) 


0.33 (0.089,0.74) 


-3.37 


1.16 


8 


13-004 


GV Tau A 


3.02(1.31,5.16) 


7.0 (6.3,7.5) 


0.55 (0.19,2.86) 


-4.11 


0.82 
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Fig. 2. Convective turnover timescales of T Tauri stars in Taurus- 
Auriga observed by XMM-Newton, plotted as a function of effec- 
tive temperature. The values have been interpolated from values 
calculated for T Ta uri stars in the Orio n Nebula Cluster (shown 
by grey squares, Preibi sch" et al.L [2005 . Preibisch, priv. comm.), 
based on stellar luminosity and effective temperature. 



4.5. Statistical analysis 

We used assumed power-law correlations between the X-ray 
activity measures and rotation pe riod, and we have u sed the 
ASURV survival analysis package (Lavalle y et al.LlT992b to per- 
form the regression. Although ASURV is d esigned to analyse 
data containing upper and/or lower limits dlsobe et all [1986) 
while our sample of T Tauri stars with measured photomet- 
ric rotation periods contained no non-detections, this enabled 
us to co mpare our results directly with those of previous stud- 
ies (e.g. IStelzer & Neuhauserl 120011: iPreibisch et all 120051) . We 
used the parametric Estimation Maximization (EM) method for 
the regression and have quoted the full range of probabilities 
of no correlation calculated by the three tests in ASURV (Cox 
Proportional Hazard, Generalized Kendall's t and Spearman's 
P)- 

We used the two-sample tests in ASURV to assess the proba- 
bilities that pairs of observed distributions of Lx/L* were drawn 
from the same distribution (one of the subsamples of stars with 
no rotational information includes an upper limit). We have 
quoted the full range of probabilities calculated by the five tests 
in ASURV (Gehan's Generalized Wilcoxon Tests using permu- 
tation variance and hypergeometric variance, Logrank Test, Peto 
and Peto Generalized Wilcoxon Test, and Peto and Prentice 
Generalized Wilcoxon Test). 



4.4. Determination of convective turnover timescales 

Convective turnover timescales m ust be determin ed empirically 
or derived from stellar models. In Preibisch et al.l d2005l) . a con- 
vective turnover timescale was derived from a full stellar model 
for each of 596 T Tauri stars in the ONC. We have determined 
the value of r conv for each T Tauri star in the XEST sample 
from the values calculated for the ONC stars (Th. Preibisch, priv. 
comm.), based on its stellar luminosity and effective temperature 
(see Fig. [2|. For multiple systems, the L* and T e ff of the pri- 
mary were used where available. The stars with T e s < 4100 K 
are essentially all fully-convect ive according to the evolution- 
ary models of ISiess et alj d2000l) and have very similar T conv of 
200-250 d. The stars in our sample with 7^ in the range 4100- 
5400 K have radiative interiors and shorter T conv , but these are 
still longer than 100 d. T con v falls steeply with further increases 
in r e ff . As the three G-type stars in the XEST sample (SU Aur, 
HD 283572 and HP Tau/G2) are much less luminous than stars 
with similar T e s in the ONC sample, all we can say is that the 
Tconv of these stars are shorter than the shortest value calculated 
for ONC stars, 11.6 d. 



5. The observed activity-rotation relations and their 
origins 

Fig. [3] (left) presents the relationship of the observed activity 
measures Lx, Lx/L+ and Fxs with rotation period for the 23 
stars in our sample with measured photometric rotation periods. 
If we consider accretors and non-accretors together as a single 
sample, as has been done in previous studies, we observe clear 
anticorrelations of Lx and Fxs with rotation period. The tests 
performed in ASURV give probabilities of less than 0.01 that 
no correlation exists. The slopes of the correlations are approxi- 
mately -1 (Table [2). However, we see no convincing anticorre- 
lation of Lx/L+ with rotation period; the tests in ASURV give 
probabilities of 0.14-0.26 that there is no correlation and the 
best-fit slope of approximately -0.4 is less than 1 .5 standard de- 
viations from zero (Table[2]i. 

Telleschi et al. I (120071) found a near-linear correlation of Lx 
with L+ in the XEST sample which they proposed was due to 
saturated activity. When they analysed the correlations of accre- 
tors and non-accretors separately, using the parametric EM al- 
gorithm in ASURV, they found similar slopes of Lx with L* for 
the two samples but that Lx was systematically 0.4 dex (a fac- 
tor 2.5) lower for accretors (see Fig . [TJ . The best-fitting relations 
were: logL x = (3O.24 + O.O6) + (1.17±O.O9)log(L*/L ) fornon- 
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Table 2. Regression parameters and their standard deviations for linear fits of observed and expected values of logLx, log(Xx/£*) 
and logfxs as a function of logP lot . A is the intercept and B is the slope. The expected activity measures were calculated for 
each star based on its stellar luminosity and whether it was an accr etor or not. The regr ession was performed using the Parametric 
Estimation Maximization method in the ASURV analysis package (lLavallev et al.L[l992h . 



Values 


logL x 


(ergs ') 






log(L x /LJ 




log F xs (erf 


;cm 1 s 


) 




A 


B 




A 


B 




A 


B 




Observed 


31.04 ±0.21 


-1.20 ± 


0.30 


-3.20 ± 


0.16 -0.37 ± 


0.24 


7.53 + 0.17 


-1.04 + 


0.25 


Expected 


30.98 ± 0.25 


-1.27 ± 


0.36 


-3.27 + 


0.10 -0.44 ± 


0.15 


7.46 + 0.17 


-1.11 ± 


0.24 



accretors andlogLx = (29.83 ±0.06) + (1.16 + 0.09) log(L*/L ) 
for accretors. 

There is the impression in Fig. [3] (left middle) that the 
accretors have consistently lower Lx/£+ at all rotation peri- 
ods, while the absence of a significant anticorrelation of Lx/£* 
with rotation period supports the concept of saturated activity. 
Nevertheless, the significant anticorrelations of Lx and Fxs with 
rotation period still suggest some kind of X-ray activity-rotation 
relation. 

Let us assume that the X-ray activity of T Tauri stars in 
Taurus- Auriga is saturated in the sense that Lx, Lx/L-t. and Fxs 
have no intrinsic dependence on rotation period, but Lx is in- 
stead characteri zed solely by the depe ndences on L+ and accre- 
tion reported by Tellesc hi et alj (1200 7). We have thus calculated 
the expected value of Lx, and hence Lx/L+ and Fxs, for each 
star in our sample based on its stellar luminosity and whether 
it is an accretor or non-accretor. Then we analysed the resulting 
correlation between each activity measure and rotation period as 
we did for the observed activity measures (Fig. [3] right). We find 
that significant anticorrelations of Lx and Fxs, with probabili- 
ties of less than 0.01 that no correlation exists, result, along with 
a shallow anticorrelation of Lx/L+. The intercepts and slopes 
of these anticorrelations are consistent with those found in the 
observed data (see Table|2]i. 

The anticorrelation of Lx with rotation period can be ex- 
plained by the the dependence of Lx on L+ and the observa- 
tion that the fast rotators in Taurus-Auriga are typically more 
luminous, while the slow rotators are typically less luminous 
(Fig- III top). The shallow anticorrelation of Lx/L* with rotation 
period can be attributed to the lower Lx of accretors at any given 
L* and the observation that the fast rotators in Taurus-Auriga 
are mainly non-accretors, while the slow-rotators are mainly ac- 
cretors (Fig. 3). The anticorrelation of Fxs with rotation period 
then arises from the shallow Lx/L-t, relation due to the observa- 
tion that the fast rotators in Taurus-Auriga have generally earlier 
spectral types, hence higher T e ff (Fig. |4] bottom), than the slow- 
rotators; recall F X s = {L x /L+)o-T* s . 

A plausible physical explanation for the dependence of ro- 
tation period on accretion and stellar lu minosities an d effec - 
tive temperatures has been described by iBouvier et al.l (Q993). 
T Tauri stars spin up if they conserve angular momentum as 
they contract toward the main-sequence. Accreting stars may be 
prevented from spinning up either by losing angular momentum 
through their strong outflows or by being magnetically coupled 
to the rotation of their accretion disks, whereas non-accreting 
stars have been allowed to spin up as expected. More massive 
stars, which have higher L+ and T e g, may originate with faster 
rotation or may spin up more quickly due to their faster evo- 
lution toward the main-sequence. We note, however, that the 
lowest-mass T Tauri stars (spectral types later than M2), which 
are largely e xcluded from our st udy, also appear to be typically 
fast rotators dHerbst et aU |2006), and other star-forming regions 
have different mass and age distributions from Taurus-Auriga, 



so the dependence of rotation on and Tgg, and hence the ob- 
served correlations of Lx, Lx/L+ and Fxs cannot be expected to 
be the same in all populations of T Tauri stars. 

An intrinsic dependence of X-ray activity on rotation period 
is not required to produce the anticorrelations of Lx and Fxs on 
rotation period observed in Taurus-Auriga, and as such anticor- 
relations have not been reported in other populations such as the 
ONC, we consider there to be no convincing evidence that they 
are intrinsic properties of the magnetic activity of T Tauri stars. 



6. Origin of the lower L x /L± of accretors 

The lower Lx/L+ of accretors compared to non-accretors in 
Taurus- Auriga has been reported before (pDamiani & Micelat 
119951; iNeuhauser et all [1991 IStelzer & Neuhauserl |2001|) . and 
has usually been attributed to their slower rotation and the action 
of an activity-rotation relationship. However, the assumption we 
made in the previous section was that accretors had an intrinsi- 
cally lower Lx/L* than non-accretors at all rotation periods. 

If we analyse the correlation of Lx/L-t with rotation period 
separately for the accretors and non-accretors, these samples, 13 
accretors and 10 non-accretors, are very small and the resultant 
fits in ASURV have large uncertainties and low significances and 
can be strongly influenced by a single data point. The probability 
that no correlation exists is approximately 0.5 for both samples. 
The slope of +0.16 + 0.19 for the non-accreting sample sug- 
gests that Lx/L-t. is not anticorrelated with rotation period, but 
the slope of -0.38 ± 0.34 for the accreting sample leaves open 
the possibility that rotation could play a role in the lower Lx/L-t, 
of accretors. The exclusion of XZ Tav£| reduces the slope for the 
accretors to -0.1 1 ± 0.34. 

We have used the larger sample of T Tauri stars with mea- 
sured projected equatorial rotational velocities, vsin/, to study 
in more detail the dependences of Lx/L+ on rotation for ac- 
cretors and non-accretors. This sample contains 30 accretors 
and 17 non-accretors. We have calculated rotation periods as 
Prot/sku = InR^/v sin i. Fig. [5] shows Lx/L-t. plotted against 
■Piot/ sin/. 

We have divided both the accretors and non-accretors into 
three subsamples: fast-rotators, with P m / sin / < 6 d, slow rota- 
tors, with P m / sin i > 6 d, and stars with no measured v sin i. The 
cumulative frequency distributions of log(Lx/^+) for these sub- 
samples are compared in Fig. [6] We can immediately see that the 
distributions for the three subsamples of accretors are at lower 
log(Lx/i*) than the distributions for the three subsamples of 
non-accretors, and that the distribution for fast-rotating accretors 



5 XZ Tau has the highest Lx/L+ of the accretors, and has the lowest 
mass and stellar luminosity of any star in our sample with a measured 
rotation period. It is moreover a binary system whose secondary was 
known to be undergoing an optical outburst that influenced the X-ray 
emission of the syste m at the time of the observation used in XEST 
jGiardinoetaTll2006h . 
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Fig. 4. The stellar luminosities (top) and effective tempera- 
tures (bottom) of T Tauri stars in Taurus-Auriga with mea- 
sured photometric rotation periods. Black triangles indicate ac- 
cretors and grey circles mark non-accretors. Filled symbols de- 
note stars observed by XMM-Newton in the present survey; 
unfilled symbols denote ad ditional stars observed by ROSAT 
dStelzer & Neuhauserl l200ll) . 



is at higher log(Lx/L*) than those for the slow-rotating accre- 
tors. The mean log(Lx/L*) of each sample is given in Table[3] 
The means of the three subsamples of non-accretors are well 
within lcr of each other, while that of the fast-rotating accre- 
tors is approximately 2cr higher than those of the other accre- 
tors. For each of the three subsamples, the accretors have a mean 
log(Lx/X*) at least 3cr lower than the non-accretors. 

We have used the two-sample tests in ASURV to assess the 
probability, P same , that the distributions of Lx/£* in pairs of 
these subsamples were drawn from the same distribution. The 
distributions of Lx/L* for the subsamples of non-accretors are 
indistinguishable from one another (P same > 0.5). The distri- 
butions of Lx/L+ for the accretors which are slow rotators and 
those with no measured vsin; are also indistinguishable from 
one another (P S ame > 0.75), and not significantly different from 
that of the fast-rotating accretors (P same = 0.09-0.35 for slow- 
rotators and 0.16-0.24 for stars with no measured vsin/). We 
find no strong evidence for a dependence of Lx/L+ on rotation 
for accreting or non-accreting T Tauri stars. 




/ 



(d) 



Fig. 5. The ratio of X-ray and stellar luminosities plotted as a 
function of rotation period for T Tauri stars in Taurus-Auriga 
observed by XMM-Newton with measured v sin i. Rotation peri- 
ods have been calculated from v sin i and the stellar radius and 
the plotted values include an unknown factor sin i and are there- 
fore upper limits. Triangles denote accretors and circles mark 
non-accretors. 



However, for each of the three subsamples, the distributions 
of Lx/L-t. of accretors and non-accretors are significantly dif- 
ferent (P same = 0.01-0.03 for fast-rotators, < 0.002 for slow- 
rotators, and 0.005-0.03 for stars with no measured v sin /)• We 
therefore find good evidence for a property of accretors other 
than their slow rotation to be the reason for their lower Lx/L*. 

If we separate the samples at P I0 J sin i — 5 d instead of 6 d, 
the distributions of Lx/L* of fast-rotating accretors and non- 
accretors are less significantly different (Psame = 0.03-0.10) but 
the distributions of Lx/X+ of accreting fast and slow rotators are 
even less significantly different (P same = 0.18-0.49). There is no 
support for a dependence of LxjL* on rotation. 

The slightly higher log(Lx/X*) of fast-rotating accretors in 
our analysis may come about from an underestimation of L* for 
some stars. Lx/L* and P m J sin / are not entirely independent 
variables as P mt sin i cc and oc y/L^. If is underesti- 
mated, a data point moves to the upper left in Fig. [5] to higher 
log(Lx/i*) and lower P ro[ / sin i. The L* of accretors can be dif- 
ficult to determine due to veiling and absorption. Three of the 
four accretors with P m / sin i < 5 d and log(Lx/£*) > -3.4 have 
suspic iously high ages of > 10 Myr in the table of Gude l et al.l 
(2007a) which suggest underestimated L+. 

The accreting stars have noticeably larger error bars than the 
non-accretors. The accretors had typically lower-quality XEST 
X-ray spectra, partly due to their lower Lx, which gives a lower 
X-ray flux from the star, and partly due to higher A^h, which 
absorbs more of this X-ray flux. However, the upper limits of 
the larger error bars are produced by models with large amounts 
of very cool plasma (To * 2 MK), absorbed by very high A^h- 
As models with such cool temperatures were not found to fit the 
spectra of either accretors or non-accretors with low A^h, we con- 
sider these unlikely to be realistic models. Such models would 
anyway point to fundamental differences in the coronae of these 
accretors and the non-accretors, in temperature rather than lumi- 
nosity. 
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Fig. 6. Cumulative frequency distributions of Lx/^* for sub- 
samples of T Tauri stars in Taurus-Auriga observed by XMM- 
Newton. Grey lines show non-accretors and black lines show ac- 
cretors. Unbroken lines mark fast rotators (P m t/ sin i < 6 d), 
dashed lines mark slow rotators (P m J sinz > 6 d), and dot- 
dashed lines mark stars with no measured v sin 



We conclude that there is no demonstrable dependence of 
Lx/L+ on rotation among T Tauri stars in Taurus-Auriga. The 
activity is consistent with being saturated at log(Lx/^*) ~ 
-3.3 for accretors and * -3.7 for non-accretors. These values 
are consistent with tho se found for T Tauri stars in the ONC 
dPreibisch et al.L 120051) . The lower Lx/Li, of accretors com- 
pared to non-accretors is not caused by their slower rotation, but 
to some other property that distinguishes accretors from non- 
accretors. The most natural interpretation is that it is a direct or 
in direct effect of the accr etion process itself. We dire ct the reader 
to lTelleschi et all d2007l) and IPreibisch et alj d2005l) for discus- 
sions of mechanisms by which accretion could reduce the X-ray 
output of accreting stars. 

Additional to these discussions, Jardi ne et alj d2006l) have 
proposed that the coronae of T Tauri stars, and hence their X-ray 
emission, are limited by pressure stripping (resulting in lower- 
mass stars having lower Lx, as observed e.g. in the ONC), but 
those of accreting T Tauri stars may be further limited by the 
encroachment of the accretion disk. As a reasonable estimate 
for the location of the inner edge of the accretion disk is at the 
corotation radius, where orbiting material rotates around the star 
with the same period as the star itself rotates, Jardine et al. fur- 
ther suggested that the wide range of Lx and Lx/L+ observed 
among the accreting T Tauri stars in the ONC could result from 
this rotation dependence. We have observed neither a wide range 
in Lx/L+ among accretors nor a positive correlation of Lx or 
Lx/L* with rotation rate. However, a directly observable depen- 
dence on rotation could be swamped by other influential param- 
eters in the model (e.g. stellar mass or coronal temperature), 
and the mass-dependence of Lx and the lower Lx/L+ of ac- 
cretors in the XEST s ample reflect those observed in the ONC 
(Tell eschi et al.L 120071) . Therefore, the Jardine et al. model ap- 
pears to be a viable description of the X-ray characteristics of 
T Tauri stars that merits further observational testing. 



Table 3. A comparison of the mean log(Lx/£*) of subsamples 
of T Tauri stars in Taurus-Auriga observed by XMM-Newton. 
Stars in the fast samples have P ml / sin i < 6 d, those in the slow 
sample have P m J sin i > 6 d, and those in the none sample have 
no measured v sin i. 



Sample 


Non-accretors 




Accretors 




N <log(L x /L*)) 


N 


<log(L x /L*)> 


Fast 


9 -3.27 ± 0.06 


17 


-3.61 ± 0.09 


Slow 


8 -3.28 ± 0.08 


13 


-3.81 ±0.09 


None 


10 -3.31 ±0.13 


17 


-3.79 ±0.10 



7. Comparison with the activity of main-sequence 
stars and that of T Tauri stars in the ONC 

Fig.|7]shows Lx/i* plotted as a function of Rossby number for 
T Tauri stars in the XEST sample with measured photometric 
rotation periods, or for which rotation periods could be esti- 
mated from vsin i as in Sect. [6] Fig.[7]a lso compares t his sa mple 
with main-sequence stars studied by iPizzolato et alj d2003l) and 
T Tau ri stars in the COUP study of the ONC dPreibisch etall 
l200l . 

The K- and M-type T Tauri stars in our sample all have val- 
ues of Ro that place them deep inside the saturated regime. The 
non-accretors with measured photometric rotation periods have 
Lx/i* and Rq consistent with the main body of saturated main- 
sequence stars. The non-accretors with rotation periods derived 
from vsin / have a little lower Lx/L+ but still within the range 
of values observed for saturated main-sequence stars. We con- 
clude that the X-ray activity of non-accreting T Tauri stars in 
Taurus-Auriga is entirely consistent with that of main-sequence 
stars showing saturated activity. 

The accretors have Rossby numbers concentrated toward the 
higher end of the range of values of the non-accretors, due to 
their typically longer rotation periods, but well within the range 
of /?o that defines saturated main-sequence stars. However, as 
we saw in Sect. [6] the Lx/L+ of accretors is significantly lower 
than that of the non-accretors and a significant fraction have 
log(Lx/i*) < -4. Such an activity level would be clearly recog- 
nised as unsaturated for a late-type main-sequence star and in- 
terpreted as evidence for a less efficient magnetic dynamo due 
to slow rotation. If T Tauri stars harbour the same kind of dy- 
namo, the T conv of these stars would need to have been system- 
atically overestimated by approximately an order of magnitude 
for th is interpretation to expla i n the l ower activity of the accre- 
tors. IWuchterl & Ts charnuter] d2003l) have presented pre-main 
sequence evolutionary models in which T Tauri stars that would 
be expected to be fully convective instead have radiative cores 
and a solar-like structure due to ongoing accretion. This would 
result in shorter convective turnover timescales than we have de- 
termined here, but detailed modelling would be required to de- 
termine how much shorter. Even if this would be the case, it 
would be the effect of accretion on the T conv that would cause 
the lower X-ray activity of accretors compared to non-accretors, 
and not their slower rotation. Our results, however, point less to 
a reduced magnetic dynamo efficiency and more to the effects of 
circumstellar material on the stellar surface and atmosphere as 
an explanation of the lower X-ray em i ssion of accretors (see e.g. 
iTelleschi et all 120071: IPreibisch et al.Ll2005l) . 

Also striking in Fig.[7]are the positions of the G-type T Tauri 
stars SU Aur, HD 283572 and HP Tau/G2, for which we could 
estimate only upper limits to T conv . Although they are among the 
fastest rotators in our sample, their short T conv moves them to the 
highest Rq, and they are almost outside the expected saturated 
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regime. One might expect similar stars with longer rotation pe- 
riods, perhaps > 5 d instead of 1-2 d, to fall outside the satu- 
rated regime and to show significantly lower Lx/L*. This may 
have been o bserved among T Tauri stars with radiative zones in 
NGC 2264 dRebull et all 120061) : while many had L x /L* at sat- 
urated levels, a significant fraction also had log(Lx/L*) < -5. 
The mean log(Lx/L*) of the stars with radiative zones was sig- 
nificantly lower than that of those without. A population of stars 
with radiative zones showing such low activity is not found in 
the XEST sample. A two-sample test of stars with and without 
radiative zones finds no significant difference in the logLx/L* 
distributions of non-accretors or accretors {P s - dme = 0.15-0.44 
and 0.17-0.21, respectively). We propose that a significant pro- 
portion of the stars with radiative zones in NGC 2264 may be 
rotating slowly enough that their activity, produced in the same 
way as solar-like main-sequence stars, is no longer saturated. 

The T Tauri stars in Taurus-Auriga are concentrated to higher 
Rossby numbers than those in the ONC. This is mainly due to 
lower-mass T Tauri stars, whose rotational properties are poorly 
studied in Taurus-Auriga, but which form a large population 
in the ONC and are mostly fast-rotating. The ONC stars have 
a larger scatter in Lx/i* (Fig. Q and the ONC accretors and 
non-accretors have much more overlap in Lx/L+ (not shown 
in Fig. |7]i than those in Taurus-Auriga. This could indicate that 
younger T Tauri stars experience greater X-ray variability, or it 
could reflect difficulties in estimating the L x and L+ due to ab- 
sorption. There are a number of non-accretors in the COUP data 
with log(Lx/^*) < -4 whose Rossby numbers indicate they 
should show saturated emission by analogy with solar-like main- 
sequence stars. Accretion cannot be the cause of lower activity 
in non-accreting T Tauri stars, so such stars are potentially very 
interesting for investigating the differences between the X-ray 
activity of T Tauri stars and late-type main-sequence stars. 

8. Comparison with previous results 

The anticorrelations of Lx, Fxs and Lx/£* with rotation 
period that we have found for the combined sample of 
acc reting and non-accre ti ng st ars differ from those found 
by IStelzer & Neuhauseri d2001l) using ROSAT PSPC data. 
Stelzer & Neuhausei] d2001l) found steeper power-law slopes for 
all three activity measures, approximately -1.5 for Lx and 
Lx/L* and approximately -2 for Fxs, and a highly significant 
anticorrelation of Lx/L* with rotation period. 

We have made a careful reexamination of the data used by 
IStelzer & Neuhauseri d2001l B. Stelzer, priv. comm.) to try to 
understand the differences between the two results. This has 
revealed systematic underestimation of Lx and overestimation 
of L+, both of which were stronger for accretors than for non- 
accretors. This caused a particularly strong underestimation of 
Lx/L+, and hence Fx, for some accretors with respect to non- 
accretors. Because the accretors have generally longer rotation 
periods than the non-accretors, this led to exaggerated slopes in 
the anticorrelations of the activity measures with rotation pe- 
riod. Significant correlations of Lx and Fxs, and perhaps even 
Lx/i*, with rotation period are nonetheless to be expected as 
the tendency for faster rotators to have high er and T e g and to 
be no n-accreting is also true in the larger Stelzer & Neuhauser 
4200JJ) sample (Fig.©. 

The underestimation of Lx was primarily due to insufficient 
accounting of the absorption of X-ray flux by interstellar and 
circumstellar material in the line of sight to the star. The ROSAT 
PSPC had, in principle, the spectroscopic capability to measure 
and account for the absorbing column density, A^h- However, its 



soft energy bandpass of 0.1-2.4 keV and low sensitivity, espe- 
cially at harder energies, combined with mo dest exposure times, 
provid ed spectra of insufficient quality for Stelz er & Neuhauseri 
(2001) to perform spectral fitting for each star. They instead con- 
verted observed count-rates to X-ray fluxes using a conversion 
factor that aimed to account for absorption through a dependance 
on hardness ratio. Our calculations using the PSPC on-axis spec- 
tral response in XSPEC have found that the maximum value 
of this conversion factor would account for absorption only for 
A^h < 2.5 x 10 20 ctrT 2 , while almost all T Tauri stars in Taurus- 
Auriga are observed through greater absorbing columns. Fig. [8] 
shows that 42 of the 45 stars detected in both the XEST and 
the ROSAT PSPC study have higher L x in the XEST stud>p, and 
that the factor of underestimation in lStelzer & Neuhauseri ( 200 1 ) 
is dependent on A^h in the manner expected and can exceed an 
order of magnitude. As accreting stars are generally surrounded 
by more circumstellar material than non-accretors, they are typ- 
ically observed through higher A^h and therefore their Lx was 
typically underestimated by a larger factor. 

The overestimation of L* in IStelzer & Neuhauseri (2001) 
ca me about from the use of the bolometric luminosities listed 
bv lKenvon & Hartmannl d 19951) . as these include emission from 
circumstellar material (mostly in the infrared) as well as the out- 
put from the star itself. While the greater amount of circumstellar 
material around accretors led to a greater underestimation of Lx 
with respect to non-accretors, so it also led to a greater overesti- 
mation of their L* . 

A secondary rea s on fo r the underestimation of Lx/L* in 
IStelzer & Neuhauseri d200ll) was their treatment of multiple sys- 
tems. The observed Lx was divided equally between the com- 
ponents, whereas we now know that Lx is approximately pro- 
portional to L+ so this approach typically underestimates the Lx 
of the primary c omponent. At the same time , the bolometric lu- 
minosities from iKen von & Hartmannl d 19951) included emission 
from all components of multiple systems, and so overestimated 
the L+ of the primary. 

These difficulties demonstrate some of the challenges faced 
by studies of the X-ray activity of T Tauri stars. While the higher- 
energy bandpasses and higher sensitivities of XMM-NewtorCs 
EPIC and Chandra's ACIS detectors have made it easier to ac- 
count for absorption than was the case for the ROSAT PSPC, 
there are still cases where A^h, and hence Lx are very poorly 
constrained. Uncertainties in Lx that account for uncertainties 
in the spectral fitting are rarely shown in the literature. Could 
the surprisingly large scatter in the COUP data be attributed to 
such uncertainties? The determination of the spectral type and 
extinction of accretors is complicated by veiling, the filling in 
of photospheric absorption lines by continuum emission from 
heated accreting material, which can make also quite uncer- 
tain. The division of the observed X-ray flux among the com- 
ponents of unresolved multiple systems is always somewhat ar- 
bitrary, and even stars which are not recognised to be multiples 
may yet have undiscovered companions. There remains room for 
improvement in studies of the X-ray emission of T Tauri stars. 
While the usefulness of precision measurements is limited, be- 
cause the Lx of any individual T Tauri star is likely to be variable 
by factors of 2-3 on timescales of weeks-years and by larger fac- 
tors during flaring on shorter timescales, it is important to avoid 
systematic problems, for example between accreting and non- 
accreting stars. 



This includes also stars with spectral types later than M3. 
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N H (10 22 cm 2 ) 

Fig. 8. L ogarithm of the ratio of X-ray luminosities in the ROSAT 
study of Stelzer & Neuhauser) d2001l) and multiplicity-corrected 
X-ray luminosities in the present (XEST) study, plotted as a 
function of absorption column density, A^h- Black triangles and 
grey circles mark accreting and non-accreting T Tauri stars, re- 
spectively. Lines show the expected ratio as a function of A^h for 
three spectral models covering the range of temperature consid- 
ered in the XEST survey (log Tq = 6.3, dashed; 7.0, full; 7.5, 
dot-dashed, with f3 = -1 in each case), assuming a count-to- 
flux conversion factor of 1.0 x 10 11 erg cm *^s 1 per ROSAT 
PSP C counts' 1 was used to c alculate the ROSAT luminosi- 
ties. IStelzer &~N euhauserl d200ll) used this factor in cases where 
the hardness ratio could not be measured, for example, for the 
most-heavily absorbed stars which were not detected in the 0. 1— 
0.4 keV band, while the maximum value used was 1.36 x 10~ n . 



9. Summary and Outlook 

We have used the XMM-Newton Extended Survey of the Taurus 
Molecular Cloud (XEST) to investigate the dependence of X-ray 
activity measures, Lx, Lx/Li, and Fxs, on rotation period. 

Using only those stars with measured photometric rotation 
periods, we found power-law anticorrelations of Lx and fxs 
with rotation periods that are significant at the 99 per cent level, 
and have slopes of approximately -1. However, we found no 
significant anticorrelation of Lx/Li, with rotation period, with 
a best-fitting slope of approximately -0.3, which is consistent 
with saturation. 

Although Stelzer & Neuhauser (2001) found steeper anticor- 
relations for all th ree of these activity measure s, reexamination 
of the data used bv lStelzer & Neuh auser ( 2001 ) has revealed that 
this can be explained by a systematic underestimation of Lx and 
overestimation of L* that preferentially affected accreting stars, 
wh ich are typically slow rotators . 



iTelleschi et al. 



( 2007) have shown that the Lx of T Tauri stars 
in Taurus-Auriga has a near-linear correlation with stellar lumi- 
nosity, suggesting saturation, but that accreting stars have sys- 
tematically lower Lx than non-accretors by a factor of 2.5 at 
any given L*. The anticorrelation of Lx with rotation period can 
be explained by the observation that the fast rotators in Taurus- 
Auriga have typically higher stellar luminosities than the slow 
rotators. The shallow anticorrelation of Lx/Li, with rotation pe- 
riod can be attributed to the observation that the fast rotators in 
Taurus-Auriga are mainly non-accretors while the slow-rotators 
are mainly accretors. The anticorrelation of Fxs = (Lx/L^crT^ 



with rotation period then comes about from the shallow Lx/Li, 
relation due to the observation that the fast rotators in Taurus- 
Auriga have typically higher T e s than the slow rotators. 

Previous studies which have reported the lower Lx/Li, of ac- 
cretors in Taurus-Auriga have usually attributed it to the typi- 
cally slower rotation periods of accretors and an anticorrelation 
of Lx/Li, with rotation period. Using a larger sample of stars 
from the XEST sample with spectral types M3 or earlier, whose 
rotation periods were derived from v sin i, we find no evidence 
that slow -rotating accretors have significantly lower Lx/Li, than 
faster-rotating accretors, or that slow -rotating non-accretors have 
lower Lx/Li, than fast-rotating non-accretors. This is consis- 
tent with saturated emission. However accretors were found to 
have significantly lower Lx/Li, than non-accretors, whether they 
were fast-rotating, slow -rotating, or had no measured v sin ;. The 
lower Lx/Li, of accretors compared to non-accretors is therefore 
not due to their slower ro tation. 

iPizzolato et al.l d2003l) have shown that the Lx/Li, of solar- 
like main-sequence stars of spectral types G to early M is de- 
termined by the Rossby number, the ratio of rotation period to 
convective turnover timescale. We have estimated the convec- 
tive turnover timescale of each star in the XEST sample by 
interpolating values c alculated for T Tauri stars in the ONC 
(Preibisch et alll2005l) and hence derived the Rossby number of 
each T Tauri star in the XEST sample with rotational informa- 
tion. 

We have found that all T Tauri stars of spectral types K 
and M in the XEST sample have Rossby numbers that lie well 
within the saturated regime of main-sequence stars. The non- 
accretors have Lx/Li, entirely consistent with those of saturated 
main-sequence stars. The accretors have, however, generally 
lower Lx/Li,, and a significant fraction have log(Lx/^*) < -4. 
Although this would be interpreted as unsaturated emission if ex- 
hibited by a main-sequence star, the determined Rossby numbers 
of these accretors are approximately an order of magnitude too 
low to be in the conventional unsaturated regime. This points to- 
ward some effect of the coupling of the stellar surface and atmo- 
sphere to circumstellar material, rather than a reduced dynamo 
efficiency, as the cause of the lower X-ray output of accreti ng 
T Tauri stars (see ITelleschi et all 120071 iPreibisch et"aill2005l) . 

Identifying the process which causes this general reduction is 
an important motivation for future studies, but just one element 
in understanding the complex interaction of magnetic activity 
and the circumstellar environment. Recent studies have found, 
for instance, that dramatic increases in accretion rate can have 
quite differing effects on the observed X-ray emissi on in dif- 
ferent cases dKastner et all l2004t lAudard et all 120051) . and that 
X-ray emission ma y also arise in shoc ks, in accretion streams 
onto the star (e.g. [ Kastner et all |2002[) or in jet outflows (e.g. 
iGiidel et alll2007bl) . 

Non-accreting stars offer simpler means for understanding 
the magnetic activity and dynamo processes in pre-main se- 
quence stars. The identification of non-accretors with Lx/Li, sig- 
nificantly below the conventional saturation level and determi- 
nation of their Rossby numbers is a potentially important probe 
of their dynamo mechanism. In particular, do there exist slow- 
rotating non-accreting T Tauri stars with low activity that are 
analogous to the low-activity slow-rotating solar-like and fully- 
convective main-sequence stars and which would provide evi- 
dence for low dynamo efficiency due to slow rotation? 

Investigations of the influence of stellar or circumstellar 
properties on the X-ray emission depend on a good character- 
ization of those properties. The T Tauri stars in Taurus-Auriga 
form a relatively small but very important population for such 
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studies, because they are closeby and they and their circumstel- 
lar environments are individually well-studied. The small sample 
of just 23 stars with measured photometric rotation periods that 
we have studied here, out of a total of ~ 200 T Tauri stars in 
Taurus-Auriga, demonstrates the potential that that further opti- 
cal photometric monitoring and X-ray campaigns could unlock. 
However, there may remain still more potential in improved 
characterization of the stellar and circumstellar characteristics in 
the much larger, but more distant, sample of T Tauri stars in the 
ONC in combination with refined analysis of the already com- 
prehensive X-ray dataset obtained in the COUP campaign. It is 
also important to study pre-main sequence populations with a 
range of ages to understand the evolution of magnetic activity 
during this interesting phase of star and planet formation. 
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Fig. 3. The X-ray activity-rotation relation of T Tauri stars in Taurus-Auriga observed by XMM-Newton. The observed (left) and 
expected (right) values of X-ray activity measures are plotted as a function of rotation period: Top, X-ray luminosity; middle, its 
ratio with stellar luminosity; bottom, the surface-averaged X-ray flux. Expected values were calculated for each star based on its 
stellar luminosity and whether it is accreting or not. Black triangles indicate accretors and grey circles indicate non-accretors. Lines 
show the best-fitting power-law correlation in each case. 
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Rossby number 

Fig. 7. The ratio of X-ray and stellar luminositie s plotted as a function of Rossby number. Grey filled circles show late-type main 
sequence stars studie d by iPizzolato et alj (l2003l) : grey open squares show T Tauri stars in the Orion Nebula Cluster studied by 
Preibi sch et all d2005l) : black circles and triangles show, respectively, non-accreting and accreting T Tauri stars in Taurus-Auriga in 
the present study. The filled black symbols show stars with measured photometric rotation periods; the open black symbols have 
periods derived from v sin i measurements and are effectiv ely upper limit s to the Rossby number. The grey leftward-pointing arrows 
represent periods derived from vsin/ measurements in Pizzolatoet al. (2003). The black rightward-pointing arrows show cases 
where only upper limits to the convective turnover time could be estimated. 



